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The finite-difference time-domain (FDTD) method is employed to solve the three dimen- 
sional Maxwell equation for the situation of near-field microscopy using a sub-wavelength 
aperture. Experimental result on unexpected high spatial resolution is reproduced by our 
computer simulation. 



§1. Introduction 

Near-field scanning optical microscopy (NSOM)B'i) is a powerful tool for the 
study of nanometer features with spatial resolution of 50-100 nm. The heart of 
NSOM is a near-field probe, which is a metal-coated optical fiber tapered to sub- 
wavelength aperture. When the probe end approaches a sample surface, the object 
is illuminated and the reemitted light is collected in the near-field region of the aper- 
ture, whose diameter determines the spatial resolution of NSOM. Optical imaging 
beyond the diffraction limit is carried out by scanning the probe on the surface. In 
addition to this fundamental principle, the resolution of NSOM is also subject to the 
tapered structure of the probe. Such a behavior has been demonstrated through our 
NSOM spectroscopy of single quantum dots.crtP 

Numerical analysis of electromagnetic field in the vicinity of the aperture and 
propagation property of light in the tapered waveguide is quite advantageous for 
the understanding of experimental results. We employ the finite-difference time- 
domain (FDTD) method© in the Mur absorbing boundary condition© to solve 
the three-dimensional Maxwell equation for the same situation as the experimental 
configuration and discuss the validity of simulation results. 



§2. Calculations 



Figures 1 and 2 show the geometries of the problem. A near-field fiber probe 
with a double tapered structure collects luminescence (A = 1/im) from a quantum 
dot buried A/10 beneath the semiconductor surface. We assume the source for 
luminescence is a point-like dipole current linearly polarized along the x direction. 
The radiation caught by A/2 aperture is transported to the tapered region clad with 
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perfect conducting metal and then coupled to the ordinary waveguide (optical fiber). 
We run the simulation with time steps of cAt = A/(40\/3) until the signal intensity 
(|i? x | 2 + \E y \ 2 + |-Ez| 2 ) , which is evaluated at (0,0,3.25A), reaches steady state. 




§3. Results 

Figure 3 shows the calculated signal intensity as a function 
of the displacement of the probe from the origin. For both scans 
along x and y directions, the full width at half maximum of the 
signal (spatial resolution of NSOM) is estimated to be around 
0.25 A, which is much smaller than the aperture diameter of A/ 2. 
This performance is beyond the fundamental principle of NSOM 
and in good agreement with the experimental result. Through 
this preliminary calculation, we demonstrate that FDTD simulation is quite useful to 
understand the behavior of light in the near-field probe and to optimize its structure 
for advanced measurements. 



Fig. 1. Schematic 
picture of simula- 
tion for NSOM. 
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Fig. 2. Cross section diagram (:rz-plane at 
y = 0) of the geometry in our 3D com- 
puter simulations for the double tapered 
fiber probe. 
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3. Intensity vs. displacement from 
the origin o in Fig. [I]. Closed circles and 
squares denote the total electric field inten- 
sity along x and y direction, respectively. 
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